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The chromospheric emission of solar-type stars in the 
young open clusters IC 2391 and IC 2602 
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ABSTRACT 

In this paper we present chromospheric emission levels of the solar-type stars in the 
young open clusters IC 2391 and IC 2602. High resolution spectroscopic data were 
obtained for over 50 F, G, and K stars from these clusters over several observing 
campaigns using the University College London Echelle Spectrograph on the 3.9-m 
Anglo- Australian Telescope. Unlike older clusters, the majority (28/52) of the solar- 
type stars in the two clusters are rapid-rotators (vsini > 20 km s _1 ) with five of 
the stars being classified as ultra-rapid rotators (vsini > 100 km s _1 ). The emission 
levels in the Calcium infrared triplet lines were then used as a measure of the chro- 
mospheric activity of the stars. When plotted against Rossby number (Nr.) the star's 
chromospheric emission levels show a plateau in the emission for Log(Nft) <J -1.1 in- 
dicating chromospheric saturation similar to the coronal saturation seen in previously 
observed X-ray emission from the same stars. However, unlike the coronal emission, 
the chromospheric emission of the stars show little evidence of a reduction in emission 
(i.e. supersaturation) for the ultra-rapid rotators in the clusters. Thus we believe that 
coronal supersaturation is not the result of an overall decrease in magnetic dynamo 
efficiency for ultra-rapid rotators. 
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1 INTRODUCTION 

Solar-type stars of spectral types F, G and K arrive on the 
Zero- Age Main-Sequence (ZAMS) with a wide range of ro- 
tation rates, with many stars experiencing extremely rapid 
rotation (projected rotational velocities, usini values, of 200 
km s _1 or more). As angular momentum loss is believed 
to be dependent upon the rotation rate of a star through 
magnetic braking (e.g. iMestel fc Spruit] Il987t ). such rapid 
rotators pose something of a problem as such stars should 
lose angular momentum much faster than they appear to 
do. 

In order to account for these rapid rotators, dynamo 
theories of young solar-type stars have used dynamo satu- 
ration to slow the angular momentum loss of rapidly rotat- 
ing stars (e.g. Irwin et al.ll2007l ; iKrishnamurthi et alj[l997t 
iBarnes fc Sofia 199(f). Below a critical rotation rate (the sat- 
uration velocity) the strength of a star's magnetic dynamo 
is believed to be related to the star's rotation rate. How- 
ever, for stars rotating more rapidly than this it is theorised 
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that the strength of the star's magnetic dynamo is no longer 
dependent upon stellar rotation. As a result, the loss rate 
of angular momentum due to magnetic braking should not 
increase for solar-type stars with rotation rates greater than 
the saturation velocity. 

The main observational evidence for dynamo saturation 
comes from the coronal X-ray emission (Lx/Lboi, the star's 
X-ray luminosity divided by the star's bolometric luminos- 
ity) of young solar-type stars. When Log(Lx/Lboi) is plotted 
against a measure of rotation rate, such as vsini, it is seen to 
initially increase with vsini until reaching a value of ~ -3.0 
when a plateau in emission level is observed for stars rotat- 
ing more rapidly than vsini > 20 km s _1 . Such saturation is 
also seen in other a ctivity indicators , such as Ha and Cal- 
cium emission (e.g. ISoderblom et all 1 19931 ) and this would 
appe ar to be consistent with the effect of dyna mo satura- 
tion l|Stauffer et al.lll997l ; IPatten fc Simonlll996l) . However, 
the mechanism causing this saturation (or even if it is dy- 
namo saturation) remains unknown. 

Furthermore, for extremely rapid rotators (vsini > 100 
km s~ 1 ) the Log(Lx/Lboi) value of these stars is then seen to 
drop below the saturation level of Log(Lx/Lboi) = -3.0. This 



2 S. C. Marsden, B. D. Carter & J.-F. Donati 



effect h as been dubbed supersaturation by iProsser et al.l 
(|l996bh . 

Afthough the mechanisms causing both saturation and 
supersaturation are yet to be determined, some sugges- 
tions have been put for ward, such as centrifugai stripping 
l|jardine fc Unruhl 19991). and the complete ti tling of emit- 
ting regions ( Stepieri. Schmitt. fc Vogesll200ll ). In order to 
better constrain both saturation and supersaturation ob- 
servationally we need to determine what activity indicators 
show signs of these effects and measure the rotational veloc- 
ities at which they occur (if at all). 

The Calcium II emission of solar-type stars is an ex- 
cellent diagnostic of chromospheric emission. Traditionally 
the Ca II H and K resonance lines at 3968A and 3934A 
have been us ed as chromospheri c diagnostics and have been 
employed bv lNoves et al.l d 19841 ) to look at the relationship 
between chromospheric emission and rotation rate for solar- 
type field stars. However observations of the Ca II H and 
K lines of faint stars can be difficult, due to often low stel- 
lar flux in the blue, and the relatively poor blue response 
of many CCD d etectors. In contrast, there have been sev- 
eral studies ( e.g. Chmiclcwski 200 01; iJames fc Jeffries 1997; 



ISoderblom et al.ll 19931 ; iFoing et al.ll 19891 ; iLinskv et alJll979h 

involving the use of the Calcium II infrared triplet (IRT) 
lines at 8498A, 8542A, and 8662A, as chromospheric di- 
agnostics. These lines are strong and their location in the 
red makes them more suitable for CCD observations. The 
strongest of the IRT lines, the 8542A line, is relatively free 
from ble nds, and is largely u ncontaminated by telluric lines, 
although lChmielewskil (|200dh shows that there is some small 
contamination in the wings of the 8542A line, but not the 
core. The IRT lines have the same upper levels as the Ca II 
H and K resonance lines, but the lower levels are populated 
radiativ ely and are n ot collisionally controlled like the H and 
K lines (|Malliklll997l ). 

Chromospheric emi ssion in the 8542A line has previ- 
ously been studied by Sodcrblo m et al.l (fl9 93) (with the 
data also presented in iKrishnamurthi et al.l 1998) for young 
solar-type stars in the Pleiades cluster (age 130 ± 20 Myrs ; 
iBarrado v Navascucs, Stauffe r. fc Javawardh ana 2004]) and 
by I James fc Jeffries! (|l997l ) for the intermediate-age (~220 
Myrs) open cluster NGC 6475. Both results show chromo- 
spheric saturation occurring in the rapidly rotating stars 
with a saturation level in the 8542 A line of Log(R'g542) ~ 
-4.2, where R's542 is the flux in the 8542A line divided by 
the star's total bolometric luminosity, see § 15.41 However, 
due to the lack of ultra-rapid rotators in both these stud- 
ies, there is no way of telling if the chromospheric emission 
shows evidence of supersaturat ion akin to that of the X -ray 
data. To further the results of ISoderblom et~ai1 (l993) and 
IJames fc Jeffries! (|l997h this paper determines the chromo- 
spheric emission in the Calcium II infrared triplet lines for 
the solar-type stars in two younger clusters (IC 2391 and 
IC 2602) containing a higher proportion of rapidly rotating 
stars, including a number in the supersaturation regime. 



2 THE OPEN CLUSTERS IC 2391 AND IC 2602 

The young open clusters IC 2391 and IC 2602 are both 
southern targets located near the galactic plane. Their HIP- 



PARCOS distances of ~145 pc |van Leeuwedfl999l ) makes 

them among the closest young open clusters. 

The age of the two clusters is given by St auffer et al.l 

(|l997l ) as 30 ± 5 Myrs. This age is based on both compar- 
isons to evolutionary isochrones as well as the upper main- 
sequence turnoff age of the clust ers. More recent work by 
IBarrado v Navascues et al.l l|2004f ) has given an age for IC 
2391 of 50 ± 5 Myrs, based on the lithium depletion of 
fainter members of the cluster. This has yet to be done for 
IC 2602. 

Given the ages of these two clusters the solar-type stars 
in IC 2391 and IC 2602 should be on, or in the last phase 
of evolving to, the ZAMS. At this age, the solar- type stars 
should not have undergone any significant magnetic braking 
and thus show a large range of rotational velocities including 
some stars with ultra-rapid rotation rates (vsini > 100 km 
s _1 ) placing them within the supersaturation regime. This 
makes IC 2391 and IC 2602 ideal targets for the study of 

saturation and supersaturation. 

Photometric work on IC 2602 was done bv lBraesI (|l962l ) 
and IWhiteoakl l|l96ll ), while similar observations were un- 
dertaken for IC 23 91 b y iHogd (li"96(ih iBuscombel (|l965h . 
iPerrv fc Hill (|l969h . and lPerrv fc Bond! (|l969h among oth- 
ers. In the past decade or so there has been renewed 
interest in these two cluste r s with ROSAT X-ray obser- 
vations dRandich et al.ll"l995l ; [Patten fc Simonl IT 996) along 
with XMM-Newton observations dMarino et al" 20051 ) . as 
well as a spectroscopic survey bv lStauffer et alTl|l997l ) and 
rotation period determinatio ns for a number of IC 2602 
st ars b y Barn es et al.l 1119991) and for some IC 2391 stars 
bv lPatten fc Simonl (| 19961 ). 

The metallicities of the two clust ers has been deter- 
mined by iRandich et al.1 (|200lT ) and iD'Orazi fc Randichl 
l|2009h to be close to solar with [Fe/H] = 0.00 ± 0.01 and 
[Fe/H] = -0.01 ± 0.02 for IC 2602 and IC 2391 respectively. 

The designations of the cluster stars used throughout 
this paper are as foll o ws. F or IC 23 91, VXR: de signation 
fromlPatten fc Simonl l|l996h . H: fromlHogd dl96Ch. L: from 



iLvngal (|l96ll ). and SHJM: f rom Istauffer et al.l d 19891 ). For 
IC 2602, R: designation from Randich et al. ([1995J), B: from 
iBraesi dl962l ). and W: from IWhiteoakl dl96ll ). Many of the 
stars have multi ple designations a nd wh ere p ossible the des- 
ignati ons from IPatten fc Simonl fl996) and IRandich et al.l 
dl995l) have been used for the stars in IC 2391 and IC 2602 
respectively. 

For stars in IC 2391 the quoted p hotom e try i s 
from lPatten fc Simonl Jl996l) . ILvngal (|l96lh . iHoggl (|l960h . 
and IStauffcr e t al.l dl989l). The rotational periods are 
from IPatten fc Simonl d 19961) a nd th e previ ous vsini val- 
ues a re from Istauffer etliL | d 19891 ) and IStauffer et al.l 
d 19971). For IC 2602 the photometry is taken from 
IProsser. Randich fc Staufferl d 1996a! ) . IRandich et all (|l995l ). 
and iBarnes et al.l (119991) . The rotational periods are from 
Barn es et al.l d 19991) and the previous vsini values are from 



Stauffer et all d , 1997l ) 



Because of the similarities between the two clusters, in 
age, distance, metallicity, etc., and the relative paucity of 
stars in the clusters (compared to say the Pleiades), IC 2391 
and IC 2602 are treated as one cluster. We present the results 
from both clusters together, although stars from each cluster 
are represented by different symbols. 

Most of the information, such as names, membership, 
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positions, colours, etc., for the stars in the two clusters can 
be found in the WEBDA databastQ 



A list of the stars with 0.4 < (V-Ic)o < 1-4 observed in 
IC 2391 and IC 2602 are given in Tables[T]and[2]respectively. 



3 SELECTION OF TARGETS 

Potential members of IC 2602 were identified by 
iRandich et all (|l995l ) from ROSAT X-ray data, with a num- 
ber of active low-mass candidate membe rs of the cluster 
found . A similar approach was take n by iPatten fc Simon] 
(1 19961) for the members of IC 2391. iProsser et alj i|l996bl ) 
and lStauffer et all (|l997l ) then carried out photometric and 
spectroscopic confirmation of membership studies (along 
with vsmi measurements), and a list of active low-mass 
members of the clu sters was produced. The listing that 
IStauffer et all ll 19971) gives for IC 2391 is not as complete as 
that given bvlPatten fc Simonl l|l99fj) so we have added some 
of the lPatten fc Simon ( 19961 ) stars to the list of potential 
members. This list is still not a complete cens us of the solar- 
type m embers in the two clusters, because, as lStauffer et al.l 
1 19971 ) points out, in neither cluster does the ROSAT sur- 
vey used to determine membership cover the entire area 
of the sky over whic h cluster members are likely to occur. 
IStauffer et all jl997l ) argue that although the list is incom- 
plete it should not be overly biased (towards active stars) 
by the selection process (using X-ray emission to determine 
membership) for at least the G dwarfs in the two clusters. 
For the K dwarfs they believe that any bias should not be 
a large effect. A more comprehensive explanation of their 
selection process along with estimates of the completeness 
of the sam ple is given in Stauffer et al.1 l|l997h . 

To th elStauffer et alJ (11997!) l ist (along with additional 
stars from IPatten fc Simonl d 19961)), we have added further 
stars from the lBraesT i 19621 ) and Whiteoakl (|l96ll ) photomet- 
ric studies of IC 2602, as well as some IC 23 9 1 sta rs from 
iHoed jl960l ), lLvnga! dlQSlf). andlStauffer et~afl jl989l ), again 
not observed bv IStauffer et~ai] (jl997). The IC 2602 p hoto- 
metric studies by Braes! 1 1962h and Whiteoakl |l961) only 
went as deep as V ~ 11, meaning that stars added to the 
list from these studies were predominantly late-F/early-G 
stars. 

With the possible exception of the above mentioned ad- 
ditional stars, we believe that most of the stars in our target 
list should be cluster members. However a further member- 
ship determination was made based on radial velocity and 
Lithium line strength, see § 15.21 

For this study as many of the solar-type stars in our 
target list were observed as possible in the time available, 
without biasing the observations through colour selection 
effects. This was achieved with observations of all of the se- 
lected stars in both IC 2391 and IC 2602 with 0.4 < (V-Ic)o 
< 1.4, corresponding to ~ early-F to ~ mid-K stars. Due to 
time constraints and the wish not to bias the observations 
by only obtaining data on a few lower-mass stars, all stars 
with (V-Ic)o J> 1-4 were excluded. Thus the observations 
should constitute a relatively unbiased sample of solar-type 
stars from early-F through to mid-K. In IC 2391 32 stars 
were observed from the target list, while in IC 2602 51 stars 
were observed from the target list. 



1 WEBDA database - http://www.univie.ac.at/webda/ 



4 OBSERVATIONS 

The spectroscopic data for this project were collected over 3 
separate observing campaigns at the 3.9-m Anglo- Australian 
Telescope (AAT) located at Siding Spring Observatory in 
New South Wales, Australia. The 3 observing runs consisted 
of 3, 3, and 2 nights in March 2000, January 2001, and Febru- 
ary 2001, respectively. 

The spectroscopic data were obtained using the Uni- 
versity College London Echelle Spectrograph (UCLES). The 
detector used for the first two runs (March 2000 and January 
2001) was the Deep Depletion MITLL3 CCD with 2048 x 
4096 15 /J,m square pixels. This chip was chosen for its excel- 
lent red response and low fringing. For the third run (Febru- 
ary 2001), the chip used was the lumogen coated MITLL2A 
CCD, again with 15 /urn square pixels. This chip does not 
have the improved red response of the MITLL3, however the 
MITLL3 chip was unavailable at the time. 

Since both of the chips are larger than the unvignetted 
field of the UCLES camera, a smaller window format (2048 
x 2896 pixels) was used to reduce read out time. Using 
the 31.6 gr/mm grating, 57 orders (#55 to #111) could 
be fitted onto the detector window. Only 49 (#63 to #111) 
were reduced (the remaining orders had little or no signal 
in them) giving full wavelength coverage from ~5000A to 
~9000A. While the actual wavelength range varied slightly 
from run to run, this variation was minimal. 

With a slit- width of ~1 arcsecond the spectral reso- 
lution varied from ~44,000 to ~46,000 during the 3 runs, 
giving a velocity resolution of ~6.8 km s -1 to ~6.5 km s _1 . 

Several exposures were taken consecutively for each star 
and added together to improve the signal-to-noise and also 
in case cosmic rays affected the IRT lines. The raw frames 
were then reduced and converted into wavelength calibrated 
spectra using the ESpRIT (Echelle Spectra Reduction: an 
Intera ctive Tool) optimal extraction routines of lDonati et al.l 
d 19971 ). 



5 RESULTS 

5.1 Measuring Rotational and Radial Velocities 

The major part of this study involves investigating the re- 
lationship between chromospheric activity and the rotation 
rates of solar-type stars in IC 2391 and IC 2602. Although 
many of the stars i n the observation list have previous v sin i 
values measured bv lStauffer et al. I l|l989l,ll997h . a number of 
stars still had none. To have a consistent i; sin i determination 
for the entire observation list, it was decided to re-determine 
the values of vsini of all the observed stars. In order to do 
this the Least-Squares Deconvolution (LSD) profiles for the 
stars were ext racted from the spe ctrum of each star using 
the routines of lDonati et al.1 (119971 ). 

An LSD profile can be considered as the resultant sum 
of the several thousand photospheric lines contained within 
each Echelle sp ectrum, yielding an "average" single high 
S /N profile (see iDonati et al. I ll997l for more details). The 
resultant gain in S/N from this process was ~10 times that 
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Landscape Table 1 to go here. 
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Landscape Table 2 to go here. 
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of the S/N in the 8542A line, depending on the spectral type 
of the star observed. The LSD profiles of three low activity 
stars of differing spectral types were also produced from the 
observations. These stars were HD 16673, a Cen A, and a 
Cen B, with spectral types of F6V, G2V and K1V respec- 
tively. 

Based on the colour of the target star an inactive star 
of similar spectral type was then rotationally broadened to 
match the target star's LSD profile. Given the high S/N of 
the LSD profiles (from ~300 to ~900) the fit to the rotation- 
ally broadened profile could usually be found to within ±1 
km s -1 . This error is based on internal consistency within 
the method by observing the same star over several epochs. 
The t>sim values o f a C en A and a Cen B are given by 
IValenti fc Fischerl (|2005l ) as 2.3 km s" 1 and 0.9 km s" 1 
respectively. Given the instrumental resolution of the ob- 
servations (6.5 to 6.8 km s _1 ) we estimate that the lower 
vsini limit we could determine for our stars was ~7 km s _1 
when using a Cen A or B as a template. This is not the 
case for HD 16673. The vsini value for this star is given by 
iNordstrom et all (|2004l ) as 8 km s _1 . This means that for 
those stars using HD 16673 as a template the lower vsini 
limit was 11 km s _ . However, no slowly-rotating stars de- 
termined to be cluster members used HD 16673 as a temple 
star. The LSD profile of one of the target stars and the 
fit to a rotationally broadened low activity star is given in 
Figure QJa). A list of theusini values for the observed stars, 
determined using this LSD method (along with comparisons 
to previous measurements) is given in Tables [1] and [2] 

Using those single stars from Tables [l] and [2] that have 
previously had their tisim values determined, the relation- 
ship between the v sin i measured here using the LSD profiles 
and previous measurements (from lstauffer et al.lll989l . ll997l ) 
was determined. This relationship is shown in Figure [2] and 
using a linear fit was found to be (with 1 a errors): 



usiniLSD^ (1.04 ±0.03) x vsim pr , 



-(1.34 ±1.63). (1) 



Along with being a useful technique for the measure- 
ment of rotational broadening, the high S/N of the LSD 
profile also makes it a good way to determine if a star is a 
double- lined spectroscopic binary. Tables [1] and [2] also note 
those stars suspected of being double-lined binaries based on 
their LSD profiles. The LSD profile of one of these binaries 
is given as an example in Figure [T^b). 

All single members of a cluster are expected to have a 
similar radial velocity. During the fitting of the LSD profiles 
of the inactive stars to the target stars in IC 2391 and IC 
2602 we have also determined the radial velocity of our tar- 
get stars. These are heliocentric radial velocities which have 
then been shifted to match the centroid of the LSD profile of 
the telluric lines in each spectrum to account for any changes 
in temperature or pressure experienced by the spectrograph 
during the coarse of the observations. This should give a zero 
point error of ~ ±0.3 km s _1 for both the target and inactive 
stars. Thus given an estimated error of ~ ±1 km s _1 for the 
method we estimate that the radial velocity measurements 
are accurate to ~ ±1.5 km s _1 . The radial velocities of the 
observed stars measured through this technique are again 
given in Tables [1] and [2] The mean and standard deviation 
of the radial velocities for the two clusters (measured from 
those stars given a M membership classification, see § 15. 2|) 
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Figure 1. LSD profiles for (a) VXR66 (from IC 2391) and (b) 
R46 (from IC 2602). The profiles have been zeroed to (a) the 
radial velocity of VXR66, 15 km s — 1 and (b) the average cluster 
velocity of IC 2602, 17.4 km s — 1 . The dashed line in (a) shows the 
LSD profile of the inactive comparison star HD16673, rotationally 
broadened to match the vsini of VXR66 (vsini = 52 km s _1 ). The 
asymmetric bumps in the bottom of the LSD profile of VXR66 
are most likely due to spots on the star's photosphere. R46 is 
most likely a double-lined spectroscopic binary. 
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Figure 2. Comparison oi vsini measurements for the stars in IC 
2391 and IC 2602 (excluding those stars with«sini < 7 km s — 1 ). 
Squares represent stars from IC 2391 and diamonds IC 2602. The 
solid line represents the linear fit to the data given in Equation ffj 
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was found to be 15.9 ± 0.7 km s _1 and 17.4 ± 1.0 km s" 1 
for IC 2391 and IC 2602 respectively. 



5.2 Cluster Membership 

As mentioned in § we believe that most of the stars in 
out target list should in fact be cluster members, based on 
previous selections. However, we have made a further deter- 
mination of the cluster membership based on the strength 
of the star's Lithium 6708A line in our observations and our 
determination of the radial velocity of the star. The strength 
of the Lithium line was only determined approximately with 
respect to the strength of the nearby Ca 6717A line, and 
given as Strong, Weak, or None. A more detailed study of 
the Lithium i n the t wo clusters has been undertaken by 
IRandich et ail l| 19971 ). iRandichl l|200ll ). and iRandich et all 
( 20011 ). 

Most of the single stars in Tables Q] and [5] with strong 
Lithium appear to have radial velocities around 16 km s _1 
for IC 2391 and 17 km s" 1 for IC 2602. Thus in order to 
qualify as a member of the cluster (M classification) stars 
have to be single (i.e. no evidence of a binary in their LSD 
profile), have strong Lithium, and have a radial velocity of 16 
± 2 km s" 1 for IC 2391 stars and 17 ± 2 km s" 1 for IC 2602. 
Stars which are classified as probable members (denoted by 
a ? classification) have a strong Lithium feature and are 
either binary or if single have a radial velocity of 16 ± 9 km 
s" 1 for IC 2391 stars or 17 ± 9 km s" 1 for IC 2602 stars. All 
other stars are classified as non-members (N classification). 
This membership is given in Tables [T] and f2] In our analysis 
of the chromospheric emission of the stars (§ 15 .4f) we have 
included all single stars with either M or ? classifications. 

Since these observations there have been other determi- 
nation s of cluster membership, such as that by 
2007]) for IC 



Platais et al 



Platais et al 



2391. The membership list of 
2007) contains a number of new IC 2391 members not con- 
tained in our li st, however this was published too late for 
our study. The IPlatais et alj (|2007l ) list also has member- 
ship determinations for the stars in our list. The determi- 
nation of m embership agr e es rea sonably well with our list. 
In addition, IPlatais et alj (|2007l ) list several IC 2391 stars 
as single-lined spectroscopic binaries that we treat as single, 
however, the contribution from the secondary to the chro- 
mospheric emission should be minimal so we have included 
them. Table [3 gi ves a list of difference s between our d e signa - 
tions and that of IPlatais et alj (|2007l ). [stauffer et al.l (|1997m 
list the star R80 (IC 2602) as a possible non-member. Based 
on our measurement of its radial velocity and Lithium we 
feel it is consistent with membership, and hence we have 
included it as such. 



Table 3. List of differences in determin ation of cluste r mem ber- 
ship and binarity between our work and IPlatais et al. | J2007T) for 
the solar-type stars IC 2391. M: cluster member, N:, non-member, 
?: probable member, B: Binary. 

Star Our Platais et al. (2007) 

Name designation designation 



L32 


M 


M + 


B 


VXR02B 


N 


M 




VXR05 


? + B 


M + 


B 


VXR07 


? 


M + 


B 


VXR08 


N 


M + 


B 


VXRll 


? + B 


M + 


B 


VXR30 


N 


M + 


B 


VXR31 


? 


N 




VXR35A 


M 


? 




VXR44 


? 


M + 


B 


VXR45A 


M 


? 




VXR50A 


? 


N 




VXR52 


M 


M + 


B 


VXR62A 


? 


M 




VXR67A 


? 


M + 


B 


VXR70 


? 


M 




VXR78 


? 


N 





iPatten fc Simonl <|l996T ). The E(V-I C ) values for the two 
clusters were fo und by multiplying the E (B-V) values by 
1.25 as given in IPinsonneault et al.l l|l998l ). thus giving val- 
ues of E(V-Irj) = 0.05 and E(V-I C ) = 0.0075 (rounded to 
0.01) for IC 2602 and IC 2391 respectively. 

The luminosities of the target stars were calculated 
from: 

L , 



M bo i = BCi + I = 4.74 - 2.51og(- 



(2) 



where 4.74 is the M bo i for the Sun (|Bessell et al.lll998l ) and 
BC i is the bolometric correction in the I band calculated us- 
ing [BesseT^StringMlow| l|l993l ). The distance moduli used 
for the clusters were 5.82 and 5.81 for IC 2391 and IC 
2602 respectively ta ken from the HIPPARCOS distances of 
van Leeuwenl(|l999r i. If no (V-Ic)o value was available then: 



BCv + V 



with the BC\ 



(3) 

calculated from 



was used i nstead 
iBessell et ail l|l99Sl) . 

Masses of the observed stars were estimated from 
iD'Antona fe Mazzitellil (|l997l ) isochrones to the nearest 0.1 
solar mass (see Figure |3}. Radii were then calculated from 
Equation [2] and: 



M bol = 4.74 - 2.51og( 



T^ ff R 2 . 

A eff0 xv 



(4) 



5.3 Stellar Parameters 

The effective temperatures of the stars given in Tables [1] 
and [2] were calculated from (V-Ic)o colours and the colour 
vs. T cff relationships of lBessell. Castelli fc Plezl (|l998l ). If no 
(V-Ic)o measurement was available for the star then (B-V)o 
was used. The (V-Ic) and (B-V) values were converted to 
unreddened values by using the reddening p arameters for 
the clusters. E(B-V) = 0.04 for IC 2602 from iBraesI (|l96ll ) 
and E(B-V) = 0.006 (rounded to 0.01) for IC 2391 from 



5.4 Measuring Chromospheric Emission 

Only those stars deemed to be probable cluster members 
and probably single (or at least single-lined spectroscopic 
binaries) had their chromospheric emission measured. Sin- 
gle stars were used because of the difficulty in determining 
the relative contributions to the emission for binary compo- 
nents. Before any chromospheric analysis was carried out on 
the spectra, telluric lines in the wings of the three IRT lines 
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Figure 3. Evolutionary plots for the stars in (a) IC 2391 and 
(b) IC 2602. Only those stars which we determined as single and 
either classified as a member or a probable member (i.e. M or ? 
in Tables [T] and [2j have been included. The VXR and R has been 
re moved from the star name. Ev olutionary tracks and Isochroncs 
of iD'Antona fc Mazziteili jl997t) are shown. The solid lines are 
evolutionary tracks for 0.8 Mq to 1.5 Mq (with the thick line 
representing 1.0 Mq). The dashed lines are age isochrones for 1, 
10, 20, 30, 50, and 100 Myrs (with the thick dashed line repre- 
senting 30 Myrs). 



were first removed from both the target and inactive star 
spectra. 

As mentioned previously there have been several stud- 
ies involving the use of the IRT lines as chromospheric di- 
agnostics, however the studies have adopted slightly differ- 
ent methods of determining the chromospheric activity of 
the stars, employing various bandwidths to measure chro- 
mosphe ric emission. The meth od we have adopted is that 
used bv lSoderblom et all (|l993h in their measurement of the 
8542 A line of solar- type stars in the Pleiades, whereby the 
chromospheric emission is calculated from the emission pro- 
file of the star. 

For all three of the IRT lines the spectrum of an in- 
active star of similar spectral type was rotationally broad- 
ened to match that of the target star. The inactive star's 
spectrum was then shifted to match the radial velocity of 
the target star and was fitted to the target star's spectrum 
by matching the wings of the IRT line being examined, as 
the chromospheric emission should be limited to the core of 
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Figure 4. The emission profile for the 8542A line of R89 (IC 
2602). The dashed line is the fitted inactive star (a Cen B in this 
case) rotationally broadened to match the target star. The emis- 
sion profile (Target - Inactive) is shown shifted by -1.0 for display 
purposes. The shaded area is the measured emission equivalent 
width. 
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Figure 5. Examples of the Calcium infrared triplet lines for stars 
of decreasing activity. The stars were chosen to have similar vsini 
values and the profiles have been shifted down in 0.2 steps for 
display purposes. 



the lines. The spectrum of the inactive star was then sub- 
tracted from the target star to produce an emission profile. 
To measure the chromospheric emission from the IRT lines, 
the area under the emission profile was calculated to pro- 
duce an emission equivalent width. The process is outlined 
in Figure [4] for the 8542A line while Figure [5] gives exam- 
ples of the emission in the IRT lines for a sample of stars of 
differing activity levels. 

By subtracting the IRT line of an inactive star (taken 
with the same instrumental setup) from that of our tar- 
get stars, the photospheric contribution to the line is re- 
moved, and the resultant emission profile contains the chro- 
mospheric emission of the star. This method also removes a 
small amount of chromospheric emission, as no star is com- 
pletely inactive. The inactive stars chosen for this project 
were the same as those used to determine the vsini of the 
target stars, namely HD16673 (F6V), a Cen A (G2V), and 
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q Cen B (K1V). These stars have near solar- metallicity (HD 
16773: [Fe/H] = -0.11, a Cen A/B: [Fe/H] = 0.12 accord- 
ing to iNordstrom et all l|2004l )). Although these inactive 
stars are not zero-age main-sequence stars (it is virtually 
impossible to get an inactive zero-age main-sequence star 
for comparison) the differences in surface gravity between 
these stars and those of IC 2391/2602 should have mini- 
mal effect on the chromospheric measurements, particularly 
for rapidly-rotating stars where rotational broadening will 
dominate the line profile. 

To measure the chromospheric emission of the IRT 
lines we have determined the chromospheric emission ra- 
tio (R'irt) for each of the IRT lines. This gives the fraction 
of each star's bolometric luminosity that is emitted from 
the chromosphere in the IRT line. This is analogous to the 
chromospheric emission ratio R' hk derived f rom C a II H 
and K observations described bv lNoves et alj |l98J), where 
the prime signifies that the photospheric emission has been 
removed. R'8542, the chromospheric emis sion ratio in the 
8542A line , has been used previo usly by ISoderblom et al.l 
l| 19931 ) and ! James fc Jeffries! (|l997f l. 

To calculate R'irt for each of the lines, the emission 
equivalent widths were first converted into the surface flux 
emitted by the IRT lines (F'irt). Using the effective tem- 
peratures of the target stars from Tables [1] and [2] the surface 
flux per Angstrom was calculated for each star by in terpo- 
lating between model atmospheres from lKuruczl |l993). This 
was then multiplied by the star's Emission Equivalent Width 
to give the surface flux in each of the IRT lines (F'irt). In 
order to remove the dependence of the surface flux upon the 
star's colour, the chromospheric emission ratio in the IRT 
lines (R'irt) was then calculated for each star. This is de- 
fined as the ratio of the surface flux in the IRT line to the 
star's total bolometric emission. 

Rirt = FlRr/crT^a , (5) 

with F'irt in W m~ 2 and T cff in K. 

The chromospheric emission ratios, emission equivalent 
widths, and surface fluxes, for single members of the two 
clusters are given in Tables f3] and [S] 



5.4.I Errors in Chromospheric Emission Measures 

Possible errors in the measurement of the chromospheric 
emission come from a number of sources; errors in the mea- 
surement of the emission equivalent width, errors in de- 
termination of effective temperature, a mismatching of the 
spectral-type of the inactive stars to that of the target stars, 
and lastly the intrinsic variability of the chromospheric emis- 
sion of the star's themselves. Errors in effective temperature 
determinations of ±100 K are likely to produce errors of less 
than 10% in the measurement of the c hromospheric emission 
ratios, while ISoderblom et al.l (|l993T ) found no perceptible 
differences in the infrared triplet lines of the inactive field 
stars they studied which cover a similar colour range to our 
data, thus a mismatch of the spectral-type of the inactive 
template star should have minimal effect of the results. Thus 
the most likely sources of error to affect our results are er- 
rors in the measuring of the emission equivalent width and 
the intrinsic variability of the stars themselves. Although 
Tables [T] and [2] show that multiple observations were taken 



of each star these were mostly taken consecutively, so un- 
fortunately cannot be used to determine any measure of the 
intrinsic variability of the stars. In order to determine the 
error in the measurement of the emission equivalent width 
on the chromospheric emission ratios of these stars we have 
plotted a comparison of the chromospheric emission ratios 
for the three infrared triplet lines (see Figure [6]). We have 
then fitted simple linear fits to the data and determined the 
following relationships: 

Log(R8 542 ) = 1-07 x Log(R^ 662 ) + 0.42 (6) 

Log(R 8542 ) = 1-05 x Log(RUs) + 0.40 (7) 

Log(R^ 662 ) = 0.95 x Log(RUs) - 0.17 (8) 

Looking at the scatter around these lines of best fit 
it can be seen that the error in Log(R's542) when plotted 
against Log(R'8662), i.e. comparing the two strongest lines, 
is smaller for more active stars than for inactive stars and 
varies from <±0.1 for active stars up to ~ ±0.2 for the inac- 
tive stars in the sample. This error is actually a cumulative 
measure of the error in both of the lines plotted and thus 
the measurement error in each individual line is most likely 
less than this. However, taking into account possible errors 
in the effective temperature we believe that ±0.1 is a rea- 
sonable estimate of the error in our Log(R'iRT) values. An 
error of ±0.1 in Log(R'mT) appears to hold even for the 
ultra-rapid rotators in the clusters, however the fitting of 
the inactive star spectrum is difficult due to the rotationally 
broadened width of the chromospheric lines in these stars, 
and thus based on fitting errors we believe that the error 
in Log(R'iRT) for the ultra-rapid rotators may be slightly 
higher than this at ±0.15. 

5.5 Activity versus Rotation 

As mentioned previously the coronal X-ray emission of 
young active solar-type stars show a saturation of activity 
for rapid rotation, with the emission showing a further su- 
pcrsaturation effect for ultra-rapid rotation. Figure [7] plots 
the Log of the coronal X-ray emission (divided by the star's 
bolometric luminosity) and the Log of the chromospheric 
emission ratio for the strongest of the IRT lines, 8542A, 
against vsmi for the observed single stars in IC 2391 and 
IC 2602. The c oronal X-ray data for IC 239 1 and IC 2602 
are taken from iPatten fc Simon! ljl996h and IStauffer et al.l 
l|l997h and were not taken simultaneously with the chorom- 
spheric results presented here. This means that variabil- 
ity and flaring may influence the scatter in these plots 
(and those in Figures [9] and lll[) . The vsmi values used 
in Figure [7] are those measured in this stu dy. The coro- 
nal d ata fo r the Pleiade s were taken from IStauffer et alj 
(|l994h and iMicela et alj (I1999T) while t he ch romospheric 
data were taken from ISoderblom et al. I 1 19931 ). For NGC 
6475 both the coronal and c hromospheric data were taken 
from I James fc Jeffries! (|l997h . 

Note: For the three stars in Tables [1] and [2] that have 
values of vsmi < 7 km s _1 and are plotted in this section 
(VXR76A, Rl, and R15A), their vsini values have been as- 
sumed to be 7 km s _1 for plotting purposes and the calcula- 
tion of Rossby number if the rotational period is unknown. 

The activity levels of the target stars have also been 
determined for the two weaker IRT lines (8498A and 8662A). 
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Figure 6. Comparison of the chromospheric emission ra- 
tios for the three infrared triplet lines, (a) Log(R's542) vs 
Log(R' g662 ), (b) Log(R' 85 4 2 ) vsLog(R' 8 498), and (c) Log(R' 86 62) 
vs Log(R'g498). As in Figure [2] squares represent stars from IC 
2391 and diamonds IC 2602. The lines are simple linear fits to 
the data and are given in Equations \E\ and \E\ 



The Log of the chromospheric emission ratio versus i; sin i for 
these two lines is given in Figure [8] 

Both the chromospheric (in all three of the IRT lines) 
and the coronal emission show a saturation for rapid rota- 
tors. The coronal emission shows supersaturation for ultra- 
rapid rotators. However, the chromospheric emission shows 
little (if any) decrease for the ultra-rapid rotators (although 
the 8662A line appears to show a possible slight decrease). 

In order to remove both the uncertainty in the incli- 
nation angle as well as any mass dependency in the ve- 
locity at which saturation starts (the saturation velocity), 
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Figure 7. (a) The Log of the coronal X-ray emission divided by 
the star's bolometric luminosity and (b) the Log of the chro- 
mospheric emission ratio in the 8542A line, both versus vsini 
for the observed single members of IC 2391 and IC 2602. Again 
squares represent stars from IC 2391 and diamonds IC 2602. Filled 
and open symbols represent stars with (V-Ic)o > 0.6 and (V- 
Ic)o < 0.6 respectively. Included in the plot are results from the 
Pleiades (crosses) and NGC 6475 (pluses). For (b) all stars with 
Log(R'8542) < -5.0 are not plotted (all such stars have (V-Ic)o 
< 0.6). 



iNoves et al.l (|l984l ') plotted activity against a value called 
the Rossby number (Nr) as a measure of a star's rotation 
rate. Rossby number is defined as the star's rotational pe- 
riod divided b y a measure of the star's convective turnover 
timescale (r r ). INoves et all (|l984l ) determined r c based on 
the (B-V) colour of the target star. 

A number of the stars in Tables [1] and [2] have period 
determinations. The Ros sby number f o r thes e stars has been 
calculated using r c from lNoves et al" (1984) Equation 4: 



log(rc) = 



1.362 
1.362 



0.166x + 0.025x 5 
■ 0.14a; 



5.323a; 3 



x > 
x < 



(9) 



where x = 1 — (B-V)o and r c is in days. There 
have been later calculat ions of the convective turnover 
timescale, such as those bv | Kim fc Demarquel (Il996l ) used by 
iKrishnamurthi et al. I (|l99&jr which are more suited to pre- 
main sequence stars. The Kim fc Demarquel (| 19961) values 
for t c do vary from those ca lculated using the 



Noves et al.l 

1984) formula above. If the iKim fc Demarquel (| 1996^ ) val- 
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Figure 8. The chromospheric emission ratio in (a) the 8498A line 
and (b) the 8662A line versus vsini for our IC 2391/2602 stars. 
As with Figure[7]all stars with log(R'g498) and log(R'g662) < -5.2 
and -5.1 respectively, are not shown. Again all these stars have 
(V-Ic)o < 0.6. The symbols are the same as in Figure[7] 



ues are used, Figure [9] remains fundamentally unchanged 
except that the Rossby numbers are all shifted to lower val- 
ues and the saturation occurs at Log(NR) ~ -1.5 rather than 
the Los (Nr) ~ -1.0 s hown in Figure [9] We have decided to 
use the iNoves et al.l l| 19841 ) values as many of the previous 
observation s of saturation and su p ersaturation have used 
Soderblom et all 1 19931 ; iPatten fc Simonl Il99rj ; 



them (e.g. „ 
IStauffer et al.lll997t ) 

For those stars in Tables [1] and [2] that do not have 
measured rotational p eriods, we have used the method of 
ISoderblom et all l|l993h to estimate the Rossby number from 
theusini of the star. The equation used is: 



2ttR 



(10) 



where R is the stellar radius (in km), vsini is in km s _1 , 
and t c (in seconds) is calculated from Equation|9] Assuming 
that stellar inclination angles are randomly distributed, then 
<sini> should be 0.785, so we have multiplied Equation 1101 
by this amount. We have also done this for the calculated 
Rossby numbers for the Pleiades and NGC 6475 which we 
have plotted. 

Figure [9] plots the Log of the coronal X-ray emission 
(divided by the star's bolometric luminosity) and the Log of 
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Figure 9. (a) The Log of the coronal X-ray emission divided by 
the star's bolometric luminosity, and (b) the Log of the chromo- 
spheric emission ratio in the 8542A line, versus the Log of the 
Rossby number (Nr) for the single observed members of IC 2391 
and IC 2602. Again results from the Pleiades (crosses) and NGC 
6475 (pluses) are included and squares represent stars from IC 
2391 and diamonds IC 2602. Filled and open symbols represent 
stars with known periods and stars whose Rossby number had 
to be estimated from Equation 1101 respectively. The black line 
is just an estimate of the saturation level for Log(NpJ < -1.1, 
while for Log(Na) > -1.1 the line is a simple linear fit to the IC 
2391/2602 data points with a fixed endpoint of Log(Nri) = -1.1 
and Log(R' 8 542) = -4.2 



the chromospheric emission ratio for 8542A line against the 
Log of the Rossby number for our stars in IC 2391/2602 with 
Figure [10] showing how the chromospheric emission changes 
with colour (unsurprisingly showing that redder stars are 
more active). Again, in both Figures, the coronal X -ray 
data for IC239 1/260 2 are from lPatten fc Simonl (|l996t ) and 
IStauffer et al.1 (|l997T ). t he coronal and chromospheric data 
for the Pleiades are from IStauffer et all ljl994h . lMicela et al.l 
l| 19991 ). a ndlSoderblom et all (Il993l ). and the NGC 6475 data 
are from I James fc Jeffries! l| 19971 ). 

As was shown in Figures [7] and [8] both the coronal 
and the chromospheric emission show a saturation level. 
Again, the coronal emission shows evidence of supersatu- 
ration (with a decrease in activity for Log(NR.) < -1.7), 
whereas the chromospheric emission shows little evidence 
of a decrease. The chromospheric saturation level is around 
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Figure 10. The Log of the chromospheric emission ratio in the 
8542A line, versus the Log of the Rossby number (Nr) for the 
single observed members of IC 2391 and IC 2602, separated into 
colour bins. Only those stars with (V-Ic)o > 0.6 have been in- 
cluded. The symbols are the same as in Figure [9] however the 
scale of the plots is different. 



Log(R'8542) ~ -4.2. This is similar to that shown in the Ca 
II 8542A emission for the Pleiades stars by Soderblo m et al.l 
l| 19931 ) and the NGC 6475 stars by I James fc Jeffries! <|l997F 7 
The scatter that the unknown inclination angle and esti- 
mated radius adds to Figure [9] for stars without a measured 
rotation period does not appear to be significant compared 
with the inherent scatter in the emission of the stars. 

The onset of saturation appears to occur at around 
Log(NR,) ~ -1.0 for both the coronal and chromospheric 
emission. If saturation occurs at the same Rossby num- 
ber then the relationship between the coronal and chromo- 
spheric emission should be linear approaching saturation. Ig- 
noring those stars with activity around the saturation level 
(in either coronal or chromospheric emission), i.e. includ- 
ing only those stars in the box in Figure [TJJ there appears 
to be a change in the relationship between the two activ- 
ity indicators for the stars in this study, between stars with 
Lx/Lboi < 20 x 10 -5 (Log(Lx/L bo i) < -3.7) and those stars 
above this value. This change occurs below both the coronal 
and chromospheric saturation level and could indicate that 
coronal and chromospheric emission have a slightly different 
relationship for stars with low activity than for more active 
stars, rather than indicating a change in the onset of satu- 
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Figure 11. Chromospheric emission ratio in the 8542A line ver- 
sus coronal X-ray activity. Again squares represent IC 2391 stars 
and diamonds IC 2602 stars. The dashed lines show the approx- 
imate saturation levels for the activity indicators. Open symbols 
represent stars with Log(Nj^) < -1.14 and those symbols with a 
plus sign inside represent those stars with Log(Np,) < -1.62. 



ration. For stars more active than Lx/Lboi ~ 20 x 10~ 5 the 
relationship between the coronal and chromospheric emis- 
sion (up to the saturation levels) could well be linear, but 
the scatter in the data is too large to determine if this is the 
case. 

Ignoring the supersaturated stars, all but two of the 
stars within the box in Figure [TTJ (and two outside) have 
Log(Nn ) > -1.14. This indicates that the Rossby number for 
the onset of saturation (for both coronal and chromospheric 
emission) is lower than Log(Na) = -1.0 indicated in Fig- 
ure[9]and is probably closer to Log(NR) = -1.1. In addition, 
Figure [TTJ shows that all but three stars with both a small 
Rossby number, Log(Nfj,) < -1.14, and low coronal activity, 
Lx/Lboi < 60 x 10 _J , (i.e. supersaturated) have Log(NR,) < 
-1.62, with only one star having a "normal" saturated level 
of coronal activity and a Log(NR) < -1.62. This appears 
to indicate that coronal supersaturation is occurring at a 
slightly higher Rossby number than shown in Figure a), 
possibly closer to log(NR.) = -1.6. 



5.6 Colour and Age Effects 

By th e age of the Hyades ~ 625 ± 50 Myrs l|Perrvman et al.1 
1998), the rotation rate of a solar-type star can be reason- 
ably determined based on its c olour, with more m assive stars 
showing higher rotation rates (R adick et al.|[l987l ). However, 
due to the minimal amount of magnetic braking they have 
undergone, solar-type stars of the age of IC 2391 and IC 2602 
should show a wide range of rotation rates for all colours. 
Figure [12] plots v sin i against colour for the single stars in IC 
2391 and IC 2602. As expected there is a wide range oft;sini 
values (up to usini = 100 km s _1 ) in the clusters regardless 
of the star's colour, however the ultra-rapid rotators in the 
clusters are all redder than (VTc)o ~ 0.85. There is a small 
dip in the data around (VTc)o ~ 0.8, showing that there 
is a paucity of rapid rotators observed for this colour. We 
suggest that this is most likely just a statistical effect. 
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Figure 12. Measured vsini values versus colour. The symbols are 
the same as in Figure [2] 



As stars age, their activity level falls as their rotation 
rate decreases through the loss of angular momen t um vi a 
magnetic braking, i.e. iBouvier. Forestini fc Allaml (|l997t ). 
Therefore it is expected that the overall activity levels of 
clusters should progressively shift to lower values for clus- 
ters of increasing age. There are two other young clusters for 
which the chromospheric emis sion ratio in the 8542A l ine has 
been measured, the Pleiades (jSoderblom et af] [l993) at an 
age of 130 ± 20 M yrs according to the Lithium-dep letion 
boundary method dBarrado v Navascues et al.l l2004f ) and 
NGC 6475 l|james fc Jeffrieslll997f ) at an age of ~220 Myrs. 
Figure [13] gives a histogram of the chromospheric emission 
ratios for IC 2391 and IC 2602, along with the ratios for the 
Pleiades and N GC 6475 clusters. The da ta for the Pleiades 
were taken from ISoderblom et alj (1 19931) whil e the d ata for 
NGC 6475 were taken from Ijames fc Jeffries! l|l997h . Most 
stars in IC 2391 and IC 2602 are seen to have a chromo- 
spheric emission level of around R's542 ~ 2 to 5 x 10 -5 com- 
pared to the Pleiades and NGC 6475, in which the peak is 
around R's 5 42 ~ 1 to 2 x 10~ 5 . The median for IC 2391/2602 
is R'8542 ~ 4.0 x 10~ 5 , for the Pleiades the median is R'8542 
~ 2.2 x 10 -5 (although this may be skewed by the fact that 
no star has R's542 < 1.0 x 10~ 5 ) and the median is R's542 ~ 
1.7 x 10 -5 for NGC 6475. As expected the median chromo- 
spheric activity level for the clusters drops with increasing 
cluster age, with IC 2391/2602 showing a much larger spread 
in chromospheric activity compared to the older clusters. 

At the age of IC 2391 and IC 2602 the solar-type stars in 
the two clusters should have undergone no significant main- 
sequence braking and should represent the rotation rate with 
which they arrive (or will soon arrive) on the ZAMS. Fig- 
ure [14] is a histogram of the measured vsini values for the 
single stars observed in the two clusters in this study, along 
with the vsini distributions for the Pleiades and NGC 6475 . 
Again the Pleiades data are from ISoderblom et all (Il993l ) 
and the NGC 6475 data are from I James fc Jeffries! (1 19971 ). 
The number of slow rotators compared to fast rotators, us- 
ing the division ofvsini = 20 km s _1 , is approximately even 
for IC 2391/2602 with 24 stars with vsini < 20 km s _1 and 
28 stars with vsini > 20 km s _1 . Thus ~54% are rapid ro- 
tators. This is much the same as the similarly aged clus- 
ter NGC 2547 (lithium depletion boundary age = 35 Myrs, 
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Figure 13. Histogram of the distribution of chromospheric emis- 
sion ratios in the 8542A line for (a) the single members of IC 
2391 & IC 2602, (b) the Pleiades, and (c) NGC 6475. The shaded 
data in (a) represents IC 2602 while the data with no shading 
represents IC 2391. Only stars with measured vsvai have been 
included. 



I Jeffries fc Oliveirall2005h which has 50% - 60% (depending 
on the inclusion of pos sible binaries or not) of its solar -type 
stars as rapid rotators (|jeffries. Totten. fc Ja mes 2000). For 
the older clusters (Pleiades and NGC 6475) magnetic brak- 
ing should have slowed many of the rapid rotators in the 
clusters. This can be seen in their vsmi distributions, with 
the Pleiades having 67 stars with vsini < 20 km s _1 and 
only 14 stars with vsini > 20 km s" 1 , ~18% rapid rotators. 
For NGC 6475 there are 23 stars with vsini < 20 km s -1 and 
only 3 stars with vsini > 20 km s _1 , ~12% rapid rotators. 
By the age of the Pleiades the median vsini of the cluster 
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Figure 14. Histogram of the distribution of vsmi for the single 
members observed in (a) IC 2391 k, IC 2602, (b) the Pleiades, 
and (c) NGC 6475. The shading and references are the same as 
in Figure [T3l Only stars with R'8542 measures have been included. 



has decreased significantly from that of IC 2391/2602, with 
the median usini of IC 2391/2602 = 21 km s _1 , the Pleiades 
= 9 km s _1 , and NGC 6475 = 8 km s" 1 . 

When looking at the difference between the G stars (0.6 
< (B-V)o < 0.9) and the early-K stars (0.9 < (B-V)o < 1.2) 
in IC 2391/2602 the percentage of rapidly rotating G stars 
is 34% while for early-K stars it is 58%. For the Pleiades the 
percentages are 15% for the G stars and 29% for the early-K 
stars. In NGC 6475 the rapid rotators makeup 0% of the G 
stars, but 33% of the early-K stars. Thus between the ages 
of IC 2391/2602 and the Pleiades the proportion of rapidly 
rotating G and K stars drops by about half (with the propor- 
tion of rapidly rotating G stars dropping slightly more than 



for the early-K stars). By the age of NGC 6475 there are 
no rapidly rotating G stars while the proportion of rapidly 
rotating early-K stars is similar to that of the Pleiades (al- 
though this may be influenced by the small number of early- 
K stars observed in NCG 6475; 2 of the 6 observed are rapid 
rotators). 

The rotational evolution of cluster stars has been stud- 
ied by a number of people including a recent study by 
llrwin et al.l (|2007l ) which includes rotational periods for a 
number of young clusters of varying ages, including rota- 
tional da ta from IC 2391/ 2602 stars. The observations pre- 
sented in llrwin et al l (j2007T ) show that the rotational period 
of G stars falls at a similar rate to that of K stars (from a 
peak around the IC 2391/2602 age) for clusters up to ~200 
Myrs. After this age the observations show that the rotation 
rate of G stars falls quite dramatically compared to that for 
K stars. This is in good agreement with the usini results 
fro m IC 2391/2602, th e Pleiades, and NGC 6475. 

llrwin et all <l2007tl have also modelled the angular mo- 
mentum evolution of young clusters as they age, incorpo- 
rating dynamo saturation into the models. The results indi- 
cate that a mass-dependent saturation velocity is required 
to reproduce the observations (taken into account via the 
mass-dependency of the Rossby number in our results, see 
§ I5.5|l . Additionally, llrwin et al. (2007) found that while 
solid-body rotation can explain the evolution of the fastest 
rotators, core-envelope decoupling (with different rotations 
in the stellar core and envelope) appears to be required to 
explain the evolution of the slowest rotators. 



6 DISCUSSION 

6.1 Chromospheric Saturation 

Figure IHJa) gives the relationship between the coronal X- 
ray emission of the single solar-type stars in IC 2391 and 
IC 2602 and Rossby number. Figure |9jb) gives the same 
plot for the chromospheric emission ratio in the 8542A line. 
For stars with Log(Nrt) < -1.1 there is a saturation in the 
chromospheric emission at a level of arou nd Log(R's542) ~ 
-4.2. T his agrees well with t he re sults of ISoderblom et al.l 
(1993) and J amcs fc Jeffries! (|l997h , who found a saturation 
in the chromospheric emission ratio for solar-type stars in 
the Pleiades and NGC 6475 of a similar level. The saturation 
levels for the two other Calcium IRT lines at 8498A and 
8662A are shown in Figure [5] and are Log(R's498) ~ -4.4 
and Log(R'8662) ~ -4.3. This confirms that, as expected, 
there is a saturation in the chromospheric emission of these 
cluster stars similar to that observed in the X-ray data from 
this and other clusters. 

If all activity indicators show the onset of saturation 
at the same value of Rossby number, this would support 
the observed saturation being caused by the same mecha- 
nism (i.e. dynamo saturation). Figure [9] shows the onset of 
saturation is around Log(NR.) ~ -1.0 to -1.1 for both the 
coronal X-ray emission and the chromospheric. The chro- 
mospheric emission in both Ha and the 8542A line from 
the Pleiades dwarfs also show the onset of sa turation at 
around Log(N R ) ~ - 1 (ISoderblom et alj Tl993h . However, 
ICardini fc Cassatellal i|2007l) show that for the Mg II emis- 
sion in a range of stars (including those in IC 2602) the onset 
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of saturation occurs around Log(Nn) ~ -0.7 (with a satura- 
tion level of Log(FMgii) ~ 6.3). Thus the onset of saturation 
for the Mg II emission appears to occur at a slightly higher 
Rossby number than what we find f or the Ca II emission 
but m ay be within the error bars as ICardini fc Cassatellal 
(2007) calculated the R ossby numbe r empir ically themselves 
rather than using the iNoves et al.l (|l984f ) value. In other 
work, the measurement of stellar activity using AV (optical 
light curve amplitud e, used as a proxy for covering frac- 
tion of starspots) by lO'Dell et al.l (Il995l ). shows the onset 
of saturation only for Log(Na) <J -1.5, if at all (using the 

es et al. I (|l984f ) Tc values). This is long after both coronal 
and chr omospheric saturation have set in, but more recent 
work bv lKrishnamurthi et al.l (| 19981 ) has questioned the use 
of AV as an activity indicator. 

Thus we conclude that both the chromospheric and 
coronal emission appear to be saturating at a similar Rossby 
number. This supports the idea that both coronal and chro- 
mospheric saturation are caused by the same mechanism. 

What are the possi ble mechanisms behind saturation? 
iJardine fc Unruhl (l999h have suggested that the effect of 
saturation seen in X-ray emission could be the result of a 
reduction in the coronal X-ray emitting volume due to cen- 
trifugal stripping of the corona for rapid rotators. As ro- 
tation r ate increases, the c o-rota tion radius of the star de- 
creases. IJardine fc Unruhl (|l999h suggest that around the 
saturation velocity the co-rotation radius moves inside the 
X-ray emitting corona. The rise in gas pressure then breaks 
the magnetic loops to form open field regions that are dark in 
X-rays. As rotation increases the reduction in emitting vol- 
ume is balanced by a rise in the density of the corona leading 
to saturation. While this idea explains coronal (X-ray) sat- 
uration, it cannot be used to explain chro mospheric satura- 
tion. T o explain chromospheric saturation. IJardine fc Unruhl 
(1999) suggest that another mechanis m, such as th e en- 
hanced continuum emission proposed bv iDovlel (|l996l ) may 
be responsible. However, if coronal and chromospheric sat- 
uration are controlled by two different mechanisms it would 
be unlikely (but not impossible) th at they would saturate at 
the same Rossby number. Thus the lJardine fc Unruhl (| 19991 ) 
explanation would appear to be unlikely. 

Another suggest ion to explain coron al saturation has 
been put forward bv lStepien et al.l (|200lh who suggest that 
coronal saturation is the result of the coronal emitting re- 
gions of a star being completely filled. This, we assume, 
would similarly be happening to the chromospheric emis- 
sion regions to account for chromospheric saturation. In this 
scenario the onset of saturation could conceivably occur at 
similar Rossby number for both coronal and chromospheric 
emission, if both are caused by strong magnetic fields erupt- 
ing through the stellar surface and extending into the chro- 
mosphere and corona. However, as explained below, this idea 
appears to fall down when trying to explain supersaturation. 

Finally it should be noted that a different explanation 
for satur ation (but not supersaturation) has been put for- 
ward by IBarnesI (|2003l ). He argues that saturation is just 
the effect of a different dynamo in operation in rapidly- 
rotating stars, a convective rather than an interface dy- 
namo. Such an explanation is supported by the finding of 
large regions of azimuthal magnetic field on the surfaces 
of young rapidly-rotating s olar-type stars (i.e. iDonati et alJ 
2003; iMarsden et al.l 120061 ). which have been explained by 



the operation of a dynamo distr i buted through out the con- 
vection zone (|Donati et al.ll2003l ). lBarne3 (2003) claims that 
saturated stars should a ll be located on what he calls the 
convective sequence fsee IBarnesI 120031 . Figure 2). However, 
this is just saying that saturated stars are more rapidly- 
rotating than non-saturated stars (i.e. rapidly-rotating stars 
are more active) which is obviou s and thus ou r results pro- 
vide no way of determining if the lBarriesl (|2003l ) explanation 
is correct. 



6.2 Chromospheric Supersaturation? 

A further complication to the problem of saturation is the 
effect of supersaturation. The X-ray emission from the IC 
2391 and IC 2602 stars, when plotted against Rossby num- 
ber in Figure |9ja) , shows the supersaturation effect occur- 
ring around Log(NR,) ~ -1.6. This has been seen in th e X-ray 
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lesser extent, in the X-ray emission of W UM a and M dwarf 
stars [jStepieri et al.ll200ll ; |james et al.ll2000l ). This study in- 
cluded measurements of the chromospheric emission of 5 
stars in IC 2391 and IC 2602 with (V-Ic)o > 0.6 andusim 
values > 100 km s~ , including 3 stars with usini values > 
200 km s _1 . The chromospheric emission for the IC 2391 and 
IC 2602 stars, Figures Ob), El and^b), show little evidence 
of a decrease in chromospheric emission for stars with v sin i 
> 100 km s" 1 or Log(N R ) < -1.6. Thus, it would appear that 
either supersaturation does not exist for the chromospheric 
emission of ultra-rapid rotators, or it occurs at a rotation 
rate significantly greater than that for coronal emission. 

It should be noted that due to the massive broadening 
of the spectral lines with such high?; sin i values, the accuracy 
of fitting of the inactive star to the target star can have an 
impact on the level of emission measured. However, a change 
of Log(R'8542) from -4.2 to -4.5 (which would be required to 
have a similar supersaturation effect to that seen in the X- 
ray emission) is about twice what we estimate the error in 
Log(R'8542) to be (see § I5.4.1[) which we consider unlikely. 

Our Cal cium infrared t riplet results are supported by 
the results o f lRandichl (|l998l ) which show no supersaturation 
in the Ha emission of ultra-rapid rotators in the 50-Myr 
old a-Persei cluste r. In addition, the Mg II stu dy of cluster 
and field stars by ICardini fc Cassatellal l|2007h also shows 
no evidence of supersaturation, however the lowest Log(NR,) 
value they have for their stars is greater than -1.6, so they 
are not really within the supersaturation regime. 

In contrast to our results, a study of a number o f active 
stars of spectral type K2 or later by IScholz et all (|2007f) 
has shown tentative evidence for supersaturation in the 
Ha emission of two stars with vsini values around 100 km 
s _1 . However, these are true pre-main sequence objects and 
the swollen convection zone of such objects may act dif- 
ferently to that of more zero- age main-sequence stars we 
hav e studied. This i s supp orted by iFeigelson et~ai] (|2003h 
and IPreibisch et al.1 (|2005l ) who find a changed (from that 
for ZAMS stars) X-ray activity - Rossby number relation- 
ship for the pre-main sequence stars in the Orion nebula. 

Like saturation , the factors causing supersaturation are 
still open to debate. Istepieri et al.l l|200lf ) have proposed that 
supersaturation is a result of the poleward migration of ac- 
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tive coronal regions, leaving the equator void of active re- 
gions and thus decreasing the coronal filling factor of the 
star. The fact that the chromospheric emission observed in 
this paper shows little evidence of supersaturation would 
count against this model being the mechanism for supersat- 
uration, as it would be expected that a poleward movement 
of active regions would also affect the chromospheric emis- 
sion. 

Such a trend in the motion of active regions is also 
not evidenced in the Doppler imaging of spot topologies 
of young active stars. Many of the stars Doppler imaged 
show rather similar spot patterns often with polar or high 
latitude spots and some lower latitude features. Comparing 
the Doppler images of stars from IC 2391/2602 in both the 
saturation (R58, vsim — 92 km s _1 ) and supersaturation 
(VXR 45A, vsini = 235 km s _1 ) regime show similar spot 
maps (|Marsden et al. 1 12004 ) with no evidence of a further 
poleward migration of active regions for the supersaturated 
star. 

The coronal stripping model of lJardine fc Unruhl (|l999l ) 
discussed previously could account for supersaturation. As 
the co-rotation radius continues to decrease with an increase 
in rotation rate, eventually enough of the coronal volume 
is fo rced open and the X-ray emission falls with rotation 
rate. IJardind l|2004h has used this model to ex plain both the 
modulation (observed bv lMarino et al.ll2003l ) and supersat- 
uration of the X-ray observations of the ultra-rapid rotator 
VXR45A. If the coronal stripping model is correct then it 
would be unlikely that chromospheric emission would show 
any evidence of supersaturation as centrifugal stripping will 
only affect the coronae of these stars. Thus the observations 
presented in this paper, showing no chromospheric supersat- 
uration, are consi s tent w ith the coronal stripping model of 
IJardine fc Unru h (1999). Unfortunately, as mentioned pre- 
viously, the fact that both the chromospheric and coronal 
emission show the onset of saturation at a similar Rossby 
number would appear to count against this model. 

Thus the mechanisms causing both saturation and su- 
persaturation remain elusive. 



7 CONCLUSIONS 

In this paper we have reported on the chromospheric emis- 
sion levels of over 50 solar- type F, G, and K stars in the 
young open clusters IC 2391 and IC 2602. Due to their 
youth, IC 2391 and IC 2602 contain a number of star that 
can be classified as ultra-rapid rotators (vs'mi > 100 km s _1 ) 
and thus are excellent targets to determine if chromospheric 
emission shows signs of supersaturation. 

We have redetermined vsini values for all our targets 
in both clusters and determined vsini values for a number 
of targets (especially those ultra-rapid rotators) for the first 
time. The clusters show a wide range of vsini values for most 
colours studied (0.4 < (V-Ic)o < 1-4) although there is a 
deficit of rapid rotators around (V-Ic)o ~ 0.8, which we 
believe has just occurred by chance. 

We have shown that chromospheric emission saturates 
for stars with Log(NR.) <C -1.1. The saturation levels found 
are Log(R' 8 542) ~ -4.2, Log(R' 8 49s) ~ -4.4, and Log(R' 86 62) 
~ -4.3, for the three Calcium IRT lines at 8542A, 8498A, 
and 8662A respectively. The chromospheric saturation level 



in the 8542A line agrees well with that for two older clus- 
ters, NGC 6475 and the Pleiades (jjames fc Jeffries! Il997l ; 
ISoderblom et al1ll993T ). Both activity indicators appear to 
be saturating at a similar Rossby number which lends weight 
to the saturation seen being caused by a single mechanism, 
which remains unknown. 

Significantly, chromospheric emission from these stars 
shows little evidence of the effect of supersaturation seen in 
the coronal X-ray emission of the same stars. Thus we be- 
lieve that X-ray supersaturation is not a result of an overall 
decrease in the efficiency of the magnetic dynamo in ultra- 
rapid rotators. 
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Table 1. Physical parameters of the solar-type stars observed in IC 2391. The first four columns give the star name, date of observation, number of observations, and signal-to-noise 
(around the 8542A line) of the combined observations. Columns 5, 6, & 7 list the photometry with a : (B-V)o estimated from (V-Ic)o, and b : (V-Ic)o estimated from (B-V)o from the 
formulas of Caldwell et al. (1993). Column 8 gives the effective temperature and Column 9 is the spectral type (from Patten & Simon 1996). Column 10 is the strength of the Lithium 
line (N: No line evident, W: Weaker than the Ca 6717A line, S: Similar or Stronger than Ca 6717A). Columns 11, 12, 13, & 14 give the rotational period ( c : represents an uncertain 
value), the radial velocity from this study, the vs'mi from this study, and the vsini from previous works. Column 15 is a note on the LSD profiles, A: some Asymmetry in the LSD profile, 
M: Multiple star, E: larger Error of ±5 km s — 1 in both Radial and Rotational Velocities, and P: Possible Binary. Column 16 has our cluster membership determination (M: Member, 



N: Non-member, and ?: 


probable member), and Columns 17, 18, c 
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mass, and radius (only calculated for sing 
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Table 2. Physical parameters of the solar-type stars observed in IC 2602. The columns are the same as in Table [T] except that the spectral type comes from d : Patten (1995), c : Braes 
(1962), and f : Whiteoak (1961). 
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Table 4. Chromospheric IRT emission of the single members observed in IC 2391. Column 1 is the star name and Column 2 is the Log of the Rossby Number for the star, where s : 
means that the Rossby Number was calculated from Equation 10. Column 3 is the Log of the X-ray emission for each star taken from Patten & Simon (1996). Columns 4, 8, & 12 are 
the emission equivalent widths of the three IRT lines (8542A, 8498A, and 8662A). Columns 5, 9, & 13 are the chromospheric emission flux in each of the IRT lines. Columns 6, 10, &: 



14 are the chromospheric emission ratios for each of the IRT lines and Columns 7, 11, & 15 


arc the Logs of the chromospheric emission ratios. 
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Table 5. Chromospheric IRT emission of the single members observed in IC 2602. The Columns are the same as in Table[4] except that the X-ray emission is taken from Stauffer et al. 
(1997). 



Star 
Name 


Log(N R ) 


L°g(L x /L bol ) 


W 8542 
(mA) 


F 8542 ( xl ° 3 ) 
(Wrri ) 


^8542 
(xl0- S ) 


L °g( R 8542) 


W 8498 
(mA) 


F 8498 < xl ° 3 ) 
(Wm~ 2 ) 


^8498 
(xlO- 3 ) 


Log(R 8 49s) 


W 8Mi2 
(mA) 


F 8662 ( xl 3 ) 
(Wm~ 2 ) 


R '8662 
(xl0" B ) 


L°g(R-8662) 


Rl 


-0.6720 E 




540 


1.9 


4.2 


-4.38 


350 


1.3 


2.9 


-4.54 


420 


1.5 


3.3 


-4.48 


R7 


-0.1719S 


-4.70 


120 


1.0 


0.7 


-5.15 


70 


0.0 


0.4 


-5.40 


70 


0.6 


0.4 


-5.40 


R8A 


-0.8175 8 




320 


1.8 


2.2 


-4.66 


180 


1.0 


1.2 


-4.92 


250 


1.4 


1.7 


-4.77 


RIO 


-1.0748 g 




470 


0.9 


3.8 


-4.42 


370 


0.7 


3.0 


-4.52 


460 


0.9 


3.8 


-4.42 


R14 


-0.8475 8 


-3.40 


570 


2.2 


4.5 


-4.35 


410 


1.6 


3.3 


-4.48 


500 


1.9 


3.9 


-4.41 


R15A 


-0.7807 


-2.92 


730 


2.0 


6.0 


-4.22 


570 


1.5 


4.5 


-4.35 


540 


1.4 


4.2 


-4.38 


R29 


-1.0274 


-3.31 


650 


1.4 


5.3 


-4.28 


400 


1.0 


3.8 


-4.42 


520 


1.1 


4.2 


-4.38 


R35 


0.4063 s 


-4.09 


390 


2.0 


2.7 


-4.57 


260 


1.4 


1.9 


-4.72 


290 


1.5 


2.1 


-4.68 


R42C 


-0.7477 8 




580 


2.2 


4.4 


-4.36 


430 


1.7 


3.4 


-4.47 


450 


1.7 


3.4 


-4.47 


R43 


-1.4513 


-3.05 


940 


2.4 


7.8 


-4.11 


(i()0 


1.5 


4.9 


-4.31 


770 


1.9 


6.2 


-4.21 


R45A 


-0.60328 


-4.29 


370 


1.9 


2.7 


-4.57 


250 


1.3 


1.8 


-4.74 


270 


1.4 


2.0 


-4.70 


R52 


-1.3977 


-3.42 


790 


2.0 


6.6 


-4.18 


450 


1.1 


3.6 


-4.44 


670 


1.0 


5.3 


-4.28 


R58 


-1.2294 


-3.22 


720 


3.5 


5.2 


-4.28 


410 


2.0 


3.0 


-4.52 


(iOO 


2.9 


4.3 


-4.37 


R59 


-1.1527 


-3.13 


780 


2.4 


6.4 


-4.19 


520 


1.6 


4.3 


-4.37 


030 


1.9 


5.1 


-4.29 


R66 


-0.5386 


-3.76 


370 


1.5 


2.7 


-4.57 


250 


1.0 


1.8 


-4.74 


280 


1.1 


2.0 


-4.70 


R68 


-1.3063 


-3.03 


870 


2.2 


7.0 


-4.15 


540 


1.4 


4.5 


-4.35 


640 


1.0 


5.1 


-4.29 


R70 


-0.4453 


-4.44 


280 


1.5 


2.1 


-4.68 


210 


1.1 


1.5 


-4.82 


240 


1.2 


1.0 


-4.80 


R72 


-0.9377 


-3.01 


770 


3.7 


5.7 


-4.24 


500 


2.4 


3.7 


-4.43 


020 


2.9 


4.5 


-4.35 


R79 


-0.2684 E 


-4.38 


160 


1.2 


1.0 


-5.00 


120 


0.9 


0.7 


-5.15 


100 


0.8 


0.6 


-5.22 


R80 


-0.5173 


-3.69 


330 


1.0 


2.7 


-4.57 


220 


0.7 


1.9 


-4.72 


240 


0.7 


1.9 


-4.72 


R83 


-0.6796 


-3.55 


520 


2.4 


3.9 


-4.41 


370 


1.7 


2.8 


-4.55 


440 


2.0 


3.2 


-4.49 


R85 


-0.5466 g 


-4.74 


ISO 


1.0 


1.0 


-5.00 


80 


0.5 


0.5 


-5.30 


110 


0.7 


0.7 


-5.15 


R88A 


-2.0748 


-3.54 


740 


1.2 


5.9 


-4.23 


530 


0.8 


3.9 


-4.41 


520 


0.8 


3.9 


-4.41 


R89 


-0.7088 


-3.54 


570 


0.9 


4.4 


-4.36 


430 


0.7 


3.4 


-4.47 


450 


0.7 


3.4 


-4.47 


R92 


-0.7777 


-3.72 


430 


1.9 


3.2 


-4.49 


270 


1.2 


2.0 


-4.70 


300 


1.3 


2.2 


-4.66 


R95A 


-1.2292 


-2.96 


090 


2.2 


5.5 


-4.26 


500 


1.6 


4.0 


-4.40 


570 


1.8 


4.5 


-4.35 


R96 


-1.1313 


-3.19 


470 


0.7 


3.5 


-4.46 


400 


0.0 


3.0 


-4.52 


480 


0.7 


3.5 


-4.46 


W79 


-0.4866 


< -4.31 


300 


1.2 


2.3 


-4.64 


220 


0.9 


1.7 


-4.77 


220 


0.9 


1.7 


-4.77 



